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c/3 , Abstract: We have analyzed the observed radial- velocity curve for the X-ray binary M33 X-7 

^ ' in a Roche model. We have analyzed the dependence between the component masses and 

the degree of filling of the optical star's Roche lobe to obtain the ratio of the masses of the 
optical star and compact object. For the most probable mass of the optical star, = 70Mq, 
I the mass of the compact object is mx = 15.55 it 3.20Mq. It has been shown that black holes 

■ with masses of nix = 15M0 and even higher can form in binaries. We present characteristic 

. evolutionary tracks for binary systems passing through an evolutionary stage with properties 



similar to M33 X-7 - type objects. According to population-synthesis analyses, such binaries 
^3 ' should be present in galaxies with masses of at least 10^^ Mq. The present number of such 

0^ . systems in M33 should be of the order of unity. We have also studied the evolutionary status 

' of the X-ray binary IC 10 X-1 with a Wolf-Rayet component, which may contain a massive 

black hole. The final stages of the evolution of the M33 X-7 and IC 10 X-1 systems should 
^ . be accompanied by the radiation of gravitational waves. 



INTRODUCTION 



The X-ray source M33 X-7 (which we will call X-7) is situated in the galaxy M33 X-7 
[1]. The periodic variability of the source was first detected by Peres et al. [2l|3], whose 
gave the first estimate for the variability period Porb = 1'^.7857 and suggested that the 
X-ray source was in a stellar binary system. 

Based on X-ray data obtained by the Einstein Observatory, ROSAT, and ASCA, 
Larson and Schulman [1] refined the period of the X-ray variability in the binary, Porb = 
3'^.4531. This value was confirmed by Dubus et al. [5]: Porb = 3*^.4535 ± 0.0005. 

In 2004, the X-ray source X-7 was identified with an optical star of magnitude 18.89m 
[7]. Photometric observations in B and V [7] confirmed that X-7 was in a binary, and 
showed that the spectral type of the optical component of X-7 is B0I-O7I and its mass 
is = 25 — 35M0 [7]. The refined period of the binary (based on XMM-Newton and 
Chandra X-ray data) Porb = 3'^.45376 ± 0.00021 was obtained [7J 

Based on the Einstein, ROSAT, and XMM Newton data, Pietsch et al. [6] derived the 
new period Porb = 3^^.453014 ± 0.000020. HST observations of the OB association HS 13, 
which contains X-7, yielded the refined spectral type for the optical component 06 III, 
and a minimum mass of m„ = 2OM0 [B]. 

Based on the spectral type of the optical component, the rate of variation of the binary 
period, the absence of an X-ray pulsar, and the analysis of the X-ray spectrum and optical 
light curves of [7], Pietsch et al.[6] suggested that X-7 is a black hole (BH). 

A detailed analysis of the photometric light curves of X-7 is presented in [8]. The 
photometric period of the binary was Porb = 3*^.4530 ± 0.0014. The radius of the optical 
component is i? = 15 — 2OM0, and its effective temperature Teff = 33000 — 47000/^ [8J. 

The first dynamical estimate of the mass of the compact object in X-7 was derived 
from the observed radial-velocity curve in |^: rrix = 15.65 ± 1.45Mq. The average mass 
for most black-hole candidates in close binary systems is near rrix ~ 8Mq. The mass of 
the compact object in X-7 significantly exceeds this estimate. Note that the mass of the 
compact object was estimated for a point-mass model in [9]. In addition, the errors of 
the observed radial- velocity curve are large (about 20-30% of the half-amplitude), which 
certainly affects the accuracy of the estimated mass of the compact object. Having in mind 
the importance of a reliable mass for the compact object in X-7 for our understanding of 
the evolution of binary stars, we have analyzed the observational data of ^ in a Roche 
model pUl [TT] . We have also analyzed the evolutionary status of the X-7 system and the 
evolution of the X-ray binary IC 10 X-1, which contains a Wolf-Rayet (WR) star and a 
possible black hole with a mass of rrix ~ 23M0. 

OBSERVATIONAL DATA 

We used the spectral data of [9] obtained from August 18 to November 16, 2006 on 
the 8.2-m Gemini North Telescope at wavelengths A = 4000 — 5000 A. The radial-velocity 
curve was derived from 22 spectra using the cross-correlation method together with a 
synthetic spectrum for A = 4150 — 4300 A and A = 4521 — 4578 A. These wavelength 
intervals contain the Hell 4200 A and Hell 4541 A lines, unblended by nebular lines. 
Orosz et al. [9] used the time for the middle of the X-ray eclipse of X-7 as the zero phase, 
Tq(HJD) = 2453967.157 ± 0.048 (this is fully consistent with the zero phase obtained 
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previously in [6], and differs from it by 95.001 ± 0.014 orbital cycles). The binary orbital 
period Porb = 3*^. 453014 [9] was obtained from the X-ray light curve. Figure [H presents 
the observed radial-velocity curve. 

Since the data of [9] are not presented in numerical form, we digitized their observed 
radial-velocity curve with an accuracy for the radial velocity of the optical component 
~ 0.6 km/s, which is substantially smaller than the mean error of the observed radial 
velocities, ~ 20 — 25 km/s. Therefore, the accuracy of our digitization of the observed 
radial-velocity curve from [9] is fully suitable for our further analysis. 

ANALYSIS FOR THE OBSERVED RADIAL- VELOCITY CURVE 

The distance between the components in the X-7 binary system is comparable to the 
radius of the optical star. The degree of Roche-lobe filling by the optical component in 
X-7 is /i = 0.777 ± 0.017 [9J. The point-mass model used to interpret the radial- velocity 
curve in [9] does not allow taking into account the proximity of the components. 

In this connection, we analyzed the radial-velocity curve for the optical component 
of X-7 in a Roche model, which enables us in a first approximation to take into account 
the tidal deformation of the star and the non-uniform temperature distribution over its 
surface, due to the effects of gravitation darkening and heating of the star's surface by X- 
rays from the relativistic object. The effect of X-ray heating in this system is insignificant. 
The bolometric luminosity of the optical component for the adopted effective surface 
temperature Tg/j = 35000K is L„ = 1.3 ■ 10^^ erg/s. The average luminosity of the X-ray 
source out of eclipse is = 5 ■ 10^'' [6]. Thus, = L^/L^ — 0.04 (in spite of the 
small amount of X-ray heating of the surface of the optical star, we including this heating 
when modeling the radial- velocity curve). Therefore, the radial- velocity curve for the 
optical component is mainly perturbed by tidal-rotational deformation of the star and 
gravitational darkening of its surface. The procedure used to model the radial-velocity 
curve is described in detail by Antokhina and Cherepashchuk [TU] and Antokhina [TT]. 
Here, we briefly review the main elements of the method. 

The binary system consists of the optical star, described by the Roche model, and 
a point-like X-ray source. The tidally deformed surface of the star was subdivided into 
~ 2600 surface elements, for each of which the emergent radiation was calculated, taking 
into account gravitation darkening, the heating of the star's surface by radiation from 
its companion (the reflection effect), and limb darkening. The effect of heating of the 
star's atmosphere by X-ray radiation from the companion was included in the procedure 
frU\ [TT] by adding the emergent and incident fluxes (without taking into account radiation 
transfer in the star's atmosphere). Since the effect of X-ray heating is very small in X-7, 
this approximation is quite acceptable. 

A model radial-velocity curve was calculated based on the absorption line profile. 
The pro- file and its equivalent width for each visible surface element with local temper- 
ature Tioc and local gravitational acceleration gioc were calculated by interpolating values 
in Kurucz's tables for Balmer lines |T2]. The integrated line profile at a given orbital 
phase was calculated by summing over the visible surface of the star, taking into account 
the Doppler effect for the local profiles, previously normalized to the continuum for each 
surface element (see [IHIIII] for more detail). The calculated integrated absorption profile 
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was used to determine the radial- velocity of the star. The radial velocity at a given orbital 
phase was calculated from the average wavelength at residual intensity levels of 1/3, 1/2, 
and 2/3 of the depth of the integrated absorption profile. The limb darkening of the stellar 
disk was calculated using the square-root law I{cos9) = Io{l — x{l — cos9) + y{l — \/cos9)), 
where 6 is the angle between the direction towards the observer and the normal to the 
surface and x and y are the limb-darkening coefficients [Ti] . 



Table 1: Parameters used in modeling the radial- velocity curves for the optical component 
in the Roche model. 



Pi-) 


3.453014 


Orbital period 




var* 


Mass of the optical star 


e 


0.0 


Eccentricity 




74.6 


Orbital inclination 




0.78-1.0 


Roche-lobe filling factor for the optical component 


/ 


1.0 


Asynchronicity factor for the rotation of the optical compo- 
nent 


Tefr(K) 


35000 


Average effective temperature of the optical component 




0.25 


Gravitational darkening coefficient 


kx 


0.04 


Ratio of the X-ray luminosity of the relativistic component to 
the bolometric luminosity of the optical component L^jL^ 


A 


1.0 


Coefficient for reprocessing of the X-ray radiation 


x\ y 


-0.186; -0.683 


Limb darkening coefficient 



* - The mass of the optical star was varied in the model calculations; [i was taken to be 
0.78, 0.85, 0.90, 0.95, 1.0. 



Since the mass of the optical star is not known exactly, we treated the masses of 
both components of the binary as parameters of the model to be fit. We carried out 
an exhaustive search over the parameters and over multiple solutions for the fit. An 
exhaustive search for the mass of the compact object was made for each mass of 
the optical component from the discrete set of values 2OM0, 3OM0, 4OM0, 5OM0, 
6OM0, 7OM0, 8OM0 and for the orbital inclination i = 74°. 6. The fit of the model to 
the observations was evaluated using the criterion. We selected the significance level 
a = 5% (see [I5] for more detail). 

In [9], the observed radial velocities were obtained from fragments of spectra containing 
the Hell 4200A and Hell 454lA. Unfortunately, fitting the radial-velocity curve based 
on the Hell model line is very time-consuming. Therefore, as we indicated above, we 
constructed the theoretical radial-velocity curve using the model H^ absorption line. Test 
calculations (Figl2]) showed no qualitative difference between the model radial-velocity 
curves obtained from the Hell 4200A and H^ lines; the average quantitative difference is 
3-4 km/s, which is substantially smaller than the errors in the observed radial velocities, 
AVr ~ 20 — 30 km/s. The calculated radial- velocity curve derived for the Hell 4200A line 
was based on the model atmosphere for the optical star [TB]. Figures 3.5 present model 
profiles of the Hell 4200 A absorption line at the two orbital phases (f) = 0.0, (j) = 0.25 
and of the H^ absorption line profile for Roche-lobe filling factors of fi = 0.78, fi = 1.0. 



3 



The result of our analysis of the observed radial-velocity curve of [9] in the Roche model 
is presented in numerical form in Table [2] and graphically in FiglHl 

Figure [7] presents the discrepancy A between the observed and theoretical radial- 
velocity curves, whose minimum corresponds to the best-fit relationship between rriy and 
rrirc. It is obvious that the Roche model for X-7 cannot be rejected at the a = 0.05 
significance level. The model can be adopted, and the optimum values for the parameters 
and their confidence intervals (errors) at the 7 = 1 — a = 0.95 confidence level can be 
estimated. 



Table 2: Dependence of component masses in X-7 for i = 74°. 6. (the errors correspond 
to the 95% confidence interval for Xm statistics, where M is the number of points in the 
radial- velocity curve). 



m„(M0) 


20 


30 


40 


50 


60 


70 


80 


m,(M0) 


7.30 ± 
1.60 


9.30 ± 
2.00 


11.0 ± 
2.40 


12.65 ± 
2.70 


14.15 ± 
2.90 


15.55 ± 
3.20 


16.95 ± 
3.40 



Our analysis using our more complex and physically more adequate model confirms 
the quantitative conclusions of [9j. Due to the small degree of Roche-lobe filling by the 
optical star, the black-hole masses obtained in our study and in [H] are quite similar: 
15.55 ± 3.2OM0 and 15.65 ± 1.45M0, respectively (for the most probable mass of the 
optical star m^, = 7OM0). 

Let us consider the dependence of the amplitude of the theoretical radial-velocity curve 
on the degree of Roche- lobe filling by the optical star. Figures H] and [5] present the model 
Hell 4200 A and absorption line profiles obtained for degrees of Roche-lobe filling 
by the optical star /i = 0.78 and 1.0. The "center of gravity" for the absorption line 
obtained for = 0.78 is more shifted towards red wavelengths, compared to the position 
for complete filling of the Roche lobe; i.e., the half-amplitude of the radial- velocity curve 
of the optical star decreases with increasing degree of Roche-lobe filling /i. 

The following calculations present the quantitative difference between the half- 
amplitude of the radial-velocity curves for the point-mass and Roche models for various 
degrees of Roche-lobe filling /i. For these calculations, we took the masses of the compact 
and optical components to be ttLx = 15.55M0 and = 7OM0. The half-amplitude of the 
radial- velocity curve for these parameters of the binary was Vc = 108.65 km/s. The differ- 
ence of the half-amplitude of the radial- velocity curve presented in Table E] is essentially 
independent of the orbital inclination of the binary (at least in the interval i = 60° — 90°). 
It is obvious that the half-amplitude of the theoretical radial- velocity curve decreases with 
increasing /i, and is 6.11 km/s smaller than the value for the point-mass model for /z = 1.0. 
The mass of the relativistic object rrix in the Roche model increases with increasing /i. 

We can see from Table [3] that, up to /i = 0.90, the difference between the half- 
amplitudes of the radial- velocity curves is close to ~ 1% of the half- amplitudes themselves. 
When the Roche-lobe filling is even higher, /i = 0.95 — 1.0, the half-amplitude of the radial- 
velocity curve decreases by ~ 3 — 6% compared to the half-amplitude for the point-mass 
model. 
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Table 3: Difference between the half-amplitude of the radial-velocity curves obtained in 
the Roche and point-mass models, VRoche and V^, as a function of the degree of Roche-lobe 
filling by the optical star fi. 





0.78 


0.85 


0.90 


0.95 


1.00 


Vc — V Roche (/) 


0.80 


1.43 


2.11 


3.24 


6.11 


m,.(M0)* 


15.55 


15.60 


15.65 


15.75 


15.95 



* - Black-hole mass obtained in the Roche model for the mass of the optical 
star rriv = TOM© . The remaining model parameters are presented in Table [TJ 



The results of fitting the observed radial-velocity curve for X-7 in the Roche model are 
close to those obtained for the point-mass model, even when the Roche lobe of the optical 
star is fully filled. For example, the difference between the mass of the compact object 
obtained in the Roche model for /x = 0.78 and in the point mass model is 0.6%, while this 
difference is 2.5% for = 1.0. This closeness of the estimates for the mass of the compact 
object obtained in the Roche and point-mass models is not universal. For example, when 
a point-mass model is used to estimate the masses of X-ray pulsars in binaries with OB 
giants, the masses of the compact objects are systematically underestimated by up to 
~ 10% (see |T7l UHl [19]). In the presence of substantial X-ray heating of the optical star, 
which appreciably changes the shape of the absorption profiles, the use of a point-mass 
model yields unreahstic results [20] . 

The similarity of the estimates for the mass of the compact object obtained in the 
Roche and point mass models for the binary X-7 is due to the fairly high mass of the 
compact object and the low degree of Xray heating, = 0.04. Thus, the fit of the 
observed radial- velocity curve for the optical star in X-7 obtained for our more adequate 
model fully confirms the quantitative conclusions drawn in [9J. 

In spite of the fact that estimates of the mass of the optical component vary from 
2OM0 to 70 Mq, we have adopted my = 70Mq as the most probable value, since this is 
the value obtained from the most complete set of observational data: a spectral analysis, 
fitting of the light curve of the optical star, and the duration of the X-ray eclipse [9J. 

Based on the mass estimate rriv = 70Mq it 6.9M0 [9], we find that the mass of the 
compact object is = 14.5 — 16.5M0. As was already noted, this black-hole mass devi- 
ates substantially from the average value rrix — 8Mq for most binaries with black holes. 
We also emphasize the Roche model fit to the radial-velocity curve can be considered 
acceptable based on the Xm criterion at the 7 = 95% confidence level (where M is the 
number of observed points in the radial- velocity curve). Further, we consider evolutionary 
scenarios that can result in the formation of systems similar to X-7. 

Recently [21J, the minimum masses of the components of the binary IC 10 X-1, which 
consists of a Wolf-Rayet star and a compact object (proposed black hole) were estimated. 
The mass of the Wolf-Rayet star is = 32.7 ± 2.6M0, and the mass of the black 
hole obtained from the observed radial-velocity curve of the optical component is rrix = 
23.1 ±2.1M0. Since no exact correspondence in time has been found between the eclipses 
of the components and the transition of the radial-velocity curve derived from the Hell 
42OOA, emission line through the 7- velocity, the conclusion that the black hole in this 
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system has a high mass must considered tentative. Nevertheless, we investigated the 
possible formation of a binary with these component masses. 

EVOLUTIONARY SCENARIOS FOR X-7 AND IC 10 X-1 

The binary X-7 is considered from an evolutionary point of view in [9]. It was con- 
cluded that current evolutionary scenarios fail to explain the existence of systems similar 
to X-7. The problem encountered by some scenarios is as follows. Since the size of the 
initially more massive star exceeds the current distance between the components, the sys- 
tem should undergo a common-envelope stage, during which the components approach 
each other. However, for the formation of a black hole with amass of about 15Mq the 
star should fill its Roche lobe immediately after the completion of the helium burning in 
the core (followed by the common-envelope stage). The radius of the star should then be 
smaller than during core burning, and so the Roche lobe should not be filled. 

We analyzed possible evolutionary tracks for X-7 and IC 10 X-1 using the "Scenario 
Machine" code. The principles behind this code have been described repeatedly, and here 
we indicate only the basic parameters used in the calculations. A detailed description of 
the "Scenario Machine" can be found in [221 [23] . 

The mass functions for black holes in X-ray binaries were calculated for various evo- 
lutionary scenarios in [21] . One of the scenarios considered in [23] predicted the existence 
of black holes with masses up to SOMq. This scenario satisfied the most important cri- 
terion for realistic evolutionary scenarios: the requirement that there exist Cyg X-l-type 
systems in the Galaxy. 

A set of evolutionary scenarios A, B, C, and W are used in |24j^, which include 
various mass-loss rates, dependences of the core mass on the initial mass of the star, and 
"massradius" dependences. These scenarios are described in detail in [231 [23] • It was 
shown in [23] that scenario B must be rejected, since it cannot explain the existence of 
a Cyg X-1- type system in the Galaxy (that is not a statistical deviation). Therefore, 
we do not consider this evolutionary scenario here. Note that scenario B is based on the 
evolutionary assumptions presented in [S]. 

Let us consider the calculated evolutionary tracks for a binary system for evolutionary 
scenarios A, C, and W. The latter two scenarios have a higher massloss rate than scenario 
A. In a first rough approximation, they can be taken to correspond to higher metallicities 
of the stars. The calculations indicate that X-7 or IC 10 X-l-type systems cannot be 
obtained from scenarios C and W. We have identified a domain of parameters in scenario 
A for which the formation of X-7 or IC 10 X-l-type systems appeared to be possible. 

A more detailed description of evolutionary scenario A is given in [221 [231 [23] • and here 
we only describe the parameters used as free parameters of the problem. In our modeling 
of possible binary precursors for X-7 and IC 10 X-1, the fraction and rate of mass loss 
in the main-sequence and WolfRayet stages, as well as the fraction of the mass of the 
pre-supernova that was enclosed by the event horizon, were varied. The star's mass-loss 
rate M is very important for two reasons: it substantially affects both the semimajor axis 

^Note that scenarios A, B, C, and W represent various assumptions concerning the physical param- 
eters of stars at various stages of their evolution, rather than specifications for types of mass exchange 
corresponding to the classification introduced by Kippenhahn and Weigert |28j . 
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of the binary and the mass of the star. The mass loss rate M in the main-sequence stage 
is described by the classical formula 



M 



aL 



(1) 



where Lis the luminosity of the star, Voo the velocity of its wind at infinity, c the speed of 
light, and a a free parameter. In scenario A, the change in the mass AM during a single 
evolutionary stage does not exceed 0.1 (M — Mcore)-, where M is the mass of the star at 
the beginning of the stage and Mcore is the mass of the stellar core. We parametrized the 
mass loss at the WolfRayet stage as AM^ij = kwRMwR, where Mwr is the maximum 
mass of the star at this stage. 

The mass of the black hole Mbh formed as a result of the explosion of a pre-supernova 
with mass Mp^eSN was calculated from the formula 



where the factor kth is the fraction of the mass of the pre-supernova enclosed by the event 
horizon during the collapse. 

The factors a and kwR took the values 0.3 (the standard value used for calculations 
in scenario A) and 0.1. The weak wind roughly imitates the metallicity of X-7, which is 
smaller than the solar value [5]; the galaxy IC 10 is also poor in metals (see, for example, 
[25]). was varied from 0.1 to 1.0. 

In the course of the population synthesis, we took the X-7 system to be a binary that 
contains a black hole with mass rriBH = 14 — 17M0 and a main sequence star at the end 
of its evolution with mass m„ = 65 — 75Mq; the orbital period of the system is Porb ^ S''. 
The IC 10 X-1 system is represented in the calculations as a binary containing a black 
hole with mass rriBH = 23 — 34Mq and a WolfRayet star with mass mwR = 17 — 35Mq; 
the orbital period is Port ^ 1.5'^. A population synthesis for 10^ binaries was carried out 
for each set of parameters of the evolutionary scenario. The results suggest that X-7- type 
binaries should exist in galaxies with masses of at least lO^^M© and a star-formation rate 
specified by the Salpeter function. The number of such systems should currently be of 
the order of unity (Fig. [H]). Figures ISlfTT] present the results for the calculations. 

Figures [8] and [TO] show the number of X-7- type and IC 10 X-l-type systems expected to 
occur in a spiral galaxy with a mass of IO^^Mq and a star formation rate specified by the 
Salpeter function. The numbers in the plots indicate curves calculated with (1) a = 0.3 
and kwB. = 0.3, (2) a = 0.1 and kwR = 0.1, and (3) a = 0.1 and kwR = 0.3. Figure [S] 
shows that there exist values for k^h for which X-7-type systems can be obtained. Taking 
into account that the star-formation rate, or at least the star-formation rate per square 
area, in M33 exceeds the average for spiral galaxies of the Local Group (see, for example, 
|25j). the existence of X-7 can be explained in scenario A. It is clear from Fig. [10] that IC 
10 X-l-type systems can obtained for a broad range of values of k^h. Since the galaxy IC 
10 has active star formation, we expect that the rate of formation of massive stars in IC 
10 is to order of magnitude comparable with this rate in the Milky Way (in spite of the 
fact that the mass of IC 10 is approximately two orders of magnitude lower than that of 
the Milky Way). A burst of star formation apparently started about 10 million years ago 
in IC 10 and is still ongoing (see, for example, [26] [27])- 



Mbh — kbhMpreSN, 



(2) 
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Figure [9] presents a characteristic evolutionary scenario resulting in the formation of 
an X-7-type system. In all the scenarios in which such a system can be obtained, the 
qualitative scheme for the evolution of the binary is very similar. At the beginning of the 
evolution, this is a massive close binary, with a mass of the primary mi = 80 — I2OM0, a 
mass of the secondary m2 = 40 — 60Mq, and an initial semimajor axis a < IOORq. The 
approximate intervals for the initial parameters of the precursors of X-7-type systems 
are also shown in Figs. [12] and [131 After the end of core hydrogen burning, the more 
massive primary fills its Roche lobe. Matter starts to flow to the secondary, as a rule, 
faster than on the nuclear timescale (stage III in [23j); the final part of the mass transfer 
can sometimes occur on close to this timescale (stage Ille in [23]). When the primary 
has lost its envelope, it turns into a Wolf Rayet star, which explodes and forms a black 
hole. The weak stellar wind in the main sequence stage, and also "correct selection" 
of the mass fraction of the pre-supernova that is enclosed by the event horizon during 
the formation of the black hole, make it possible to form a black hole of the required 
mass (around 15 — 2OM0) in a sufficiently close system. Note also that an X-7-type 
system does not undergo a common-envelope stage before its formation^. Qualitatively, 
the subsequent evolution of an X-7- type binary is as follows. The secondary (initially 
less massive) star fills its Roche lobe, leading to a stage of supercritical accretion onto 
the black hole; a common-envelope then forms, and the binary ends its existence with 
the merger of its components and the formation of a Thorne-Zhytkov object. The final 
result of the evolution of an X-7-type system is a single massive black hole. When the 
black hole moves along a spiral trajectory into the center of the non-degenerate star (the 
Thorne-Zhytkov object stage), the system can be a source of gravitational waves [55] . 

Let us now consider the characteristic evolutionary track of a binary that results in the 
formation of an IC 10 X-l-type system (Fig. [TTI) . At the beginning of the evolution, the 
mass of the primary is mi = 80 — I2OM0, the mass of the secondary is m2 = 15 — 6OM0, 
and the initial semimajor axis of the binary orbit is a ^ 170 — 200/?©. The approximate 
intervals of the initial parameters of precursors of IC 10 X-l-type systems are also shown 
in Figs. [12] and [131 Note that, since the mass of a star can exceed its initial mass due to 
mass transfer, the minimum initial mass of the secondary can be lower than the estimated 
mass for the WolfRayet star in IC 10 X- 1. After the depletion of hydrogen in its core, 
the more massive primary fills its Roche lobe. Matter starts to fiow to the secondary, as 
a rule, faster than on the nuclear timescale (stage III in [23]); the final part of the mass 
transfer occurs on a timescale close to the nuclear timescale (stage Ille in [23j). When the 
primary has lost its envelope, it turns into a WolfRayet star, which explodes and forms 
a black hole. Further, the secondary finishes its evolution on the main sequence and fills 
its Roche lobe. After a phase of supercritical accretion onto the black hole during the 
Roche-lobe filling stage, a common-envelope stage begins, during which the stars approach 
each other very closely and the non-relativistic component loses its envelope. A binary 
consisting of a black hole and a WolfRayet star is formed. In the process of collapse of 
the WolfRayet star the gamma-burst may occur in such a binary. The evolution of an 

^See [53] for more detail. The "Scenario Machine" assumes that type B systems, according to 
the classification of Kippenhahn and Weigert, form a common envelope ii q < qcr — 0.3. Here, 
q = Maccretor / Mdonor is the ratio of the masses of the primary (accretor) and secondary (donor) stars. 
For the track presented in Fig. [9] g « 0.35 at the beginning of the stage of Roche-lobe filling of the 
initially more massive primary star. 
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IC 10 X-1- type binary ends with the merging of the two black holes (remnants of the 
evolution of the system components), driven by gravitational-wave radiation, resulting in 
the formation of a single massive black hole. The merger should be accompanied by a 
burst of gravitational- wave radiation. 

CONCLUSION 

We have used the radial-velocity curve for the optical star in the M33 X-7 system [9] 
and parameters of the system derived from an analysis of the light curve [9] (/x = 0.78, 
i = 74°. 6) to construct the relationship between the masses of the black hole and 
the optical star rriy (Table [21 Fig. [6]) in a Roche model for the optical star. We also 
studied the effect of tidal deformation and gravitational darkening of the optical star on 
the mass of the black hole rrix for various degrees of filling of the Roche lobe of the optical 
star /i (Table [3l). For /i = 0.78, the impact of the proximity of the components on the 
radial-velocity curve of the optical star is relatively small, justifying the analysis used by 
Orosz et al. [9] to derived the black-hole mass m^. = 15.65 ± 1.45M0 (for rriy = 70Mq). 
Our estimated mass for the black hole for m„ = 7OM0 coincides with this value within 
the errors: m^. = 15.55 ± 3.2OM0; the error corresponds to the 95% confidence interval 
for xti statistics, where M is the number of data points in the radial-velocity curve. We 
stress that, if the degree of Roche- lobe filling by the optical star were close to unity, 
with all other factors being the same, the black-hole mass would be rrix = 15.95 ± 3.2M0, 
which exceeds the value for /i = 0.78 by ~ 3%. Therefore, it is important to use a 
more realistic Roche model that takes into account the tidalrotational deformation of the 
optical star and its gravitational darkening, rather than a point mass model, when fitting 
the radial-velocity curve. 

We also considered the evolutionary status of the X-ray binaries M33 X-7 and IC 10 
X-1, and have shown that the existence of such systems, which probably contain massive 
stellar black holes, is quite possible in the galaxies M33 and IC 10 for reasonable physical 
assumptions about the initial parameters of these systems and the mass transfer in them. 
In the evolution process the IC lOX-1 system should pass through the stage, during which 
it will consist of a black hole and the WolfRayet star. At the instant of collapse of the 
WolfRayet star the gamma-burst may occur in this binary. We conclude that the evolution 
of the IC 10 X-1 system should end with a merger of two massive black holes, which should 
result in a burst of gravitational-wave radiation. In the case of the X-7 system, a burst of 
gravitational wave radiation should be observed during the formation of a ThorneZhytkov 
object [33]. These conclusions are important for observational programs to be undertaken 
with modern laser gravitational-wave telescopes, such as LIGO, VIRGO, etc. 

Note that the detection of a black hole in M33 X-7 [9] has opened the era of searches 
for black holes in other galaxies of the Local Group. This has become possible due to the 
commissioning of large, 8-10-m, new-generation telescopes. The first detections of new 
stellar-mass black holes in other galaxies are consistent with the hypothesis that there is 
a bimodal mass distribution for relativistic objects [29]- [31]. This distribution and the 
presence of a dip in the mass distribution for neutron stars and black holes in the range 
rrix = 2 — 4M0 [21] were first deduced from the masses of about 40 objects (~ 20 neutron 
stars and ~ 20 black holes). The detection of new black holes in other galaxies will make 
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it possible to substantially increase the number of mass estimates for relativistic objects, 
and thus to improve the significance of corresponding statistical conclusions. The increase 
in the number of detected black holes in X-ray binaries will also open fresh opportunities 
for verifying new theories of gravitation [311 [32] . 
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Fig. 1: Observed and theoretical radial- velocity curve for the optical component in the 
X-ray binary M33 X-7. The points show the radial velocities of the optical star taken 
from [9]. The solid curve shows the theoretical radial- velocity curve in the Roche model 
calculated for a mass of the compact object m^^ = IS.SSM©, a mass of the optical star 
rriv = 70Mq, and the orbital inclination i = 74.6°. (the remaining parameters of the 
binary are presented in Table [1]) . 
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Fig. 2: Theoretical radial-velocity curves for the optical star, obtained in the Roche 
model with the mass of the compact object = IS.SSM©, the mass of the optical star 
rriv = 70Mq, and the orbital inclination i = 74.6° (the other parameters of the binary are 
presented in Tabled]). The solid and dashed curves represent the radial- velocity curves 
calculated from the and Hell 4200A absorption lines, respectively. 
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Fig. 3: Model integrated profiles for the Hell 4200 A absorption line of the optical 
component in X-7 at orbital phase 0.0. The profiles were obtained in the Roche model 
with m^. = 15.55M0, my = 70Mq, and i = 74°. 6. The solid and dashed curves represent 
the absorption-line profiles obtained in the Roche model for degrees of filling of the Roche 
lobe fi = 0.78 and /i = 1.0, respectively (the other parameters of the binary are presented 
in Table [1]). Both model integrated profiles were convolved with a spectrograph response 
function with FWHM=2A. 
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Fig. 4: Same as Fig|3]for orbital phase 0.25. 
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Fig. 5: Same as FigJUfor the absorption line. 
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Fig. 6: Relationship between the masses of the compact object and optical component 
of X-7, obtained from the fit of the observed radial-velocity curve from ^ in the Roche 
model. The parameters of the binary are presented in Table 1. The solid curve represents 
the central values for the mass of the black hole, while the dashed curves correspond to 
the upper and lower boundaries for the errors on the central values based on the Xm 
statistics for the 7 = 95% confidence level (M is the number of observed points on the 
radial- velocity curve). 
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Fig. 7: Dependence of the discrepancy A between the observed and theoretical radial- 
velocity curve for rriv = 70Mq. The horizontal line corresponds to the critical value for 
xii a — 0.05 significance level. 
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Fig. 8: Number of M33 X-7-type systems for a spiral galaxy with mass 10^^ Mq and a star- 
formation rate specified by the Salpeter function. The numbers indicate curves calculated 
for various sets of evolutionary parameters (see the text). The vertical axis plots the 
expected absolute number of M33X-7- type systems in the Galaxy. The horizontal axis 
plots the fraction of the mass of the pre-supernova that is enclosed by the event horizon at 
the instant when the black hole is formed. In the case of curve "1" {a = 0.3, kwR = 0.3, 
ksH = 0.3), the maximum expected number of M33 X-7-type systems is 0.37, which is of 
the order of magnitude of unity. 



20 



CraflMsi M M a T 




Fig. 9: A characteristic evolutionary scenario resulting in the formation of an M33 X-7- 
type system. The notation for the evolutionary stages is (see [23]): for more detail): I main 
sequence stage, III, Ills stage of Roche lobe filling, WR Wolf-Rayet stage, BH black hole, 
SBH black hole with a supercritical accretion rate, SN supernova explosion, CE common 
envelope stage, TZ Thorne-Zhytkov object. The mass of the primary and secondary stars 
M (in M0), the semimajor axis a (in Rq), and the age of the system T (in millions of 
years) are presented for the onsets of the corresponding stages and, in the case of the 
supernova explosion, for the time just before the explosion. The free parameters of the 
problem for this scenario are a = 0.3, k^ji = 0.3, and = 0.3. 
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Fig. 10: Same as Fig. 8 for an IC 10 X-l-type system. 
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Fig. 11: A characteristic evolutionary scenario resulting in the formation of an IC 10 
X-l-type system. The notation for the evolutionary stages is (see [23]): for more detail): 
I main sequence star; II star has ended its evolution on the main sequence but does not 
fill its Roche lobe; III, Ille, and Ills stage of Roche lobe filling; WR Wolf-Rayet stage; BH 
black hole; SBH black hole with a supercritical accretion rate; SN supernova explosion; 
CE common envelope stage. The mass of the primary and secondary stars M (in Mq), the 
semimajor axis a (in Rq), and the age of the system T (in millions of years) are presented 
for the onsets of the corresponding stages and, in the case of the supernova explosion, for 
the time just before the explosion. The free parameters of the problem in this scenario 
are a = 0.1, kwR = 0.1, and k^h = 0.5. 
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Fig. 12: Initial parameters of systems that evolve into binaries of the type considered. The 
vertical axis plots the initial mass of the primary (which is initially more massive) and the 
horizontal axis the initial semimajor axis of the system. The triangles and circles indicate 
binaries that evolve to form M33 X-7-type and IC 10 X-l-type systems, respectively. The 
free parameters of the problem in this case are a = 0.1, = 0.1, k^h = 0.3. 
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Fig. 13: Initial parameters of systems that evolve into binaries of the type considered. 
The vertical axis plots the initial mass of the primary (which is initially more massive) 
and the horizontal axis the component-mass ratio q = M2/M1 < 1 at the beginning of 
the system's evolution. The notation is the same as in Fig. 12. The free parameters of 
the problem in this case are a = 0.1, kwR = 0.1, kbh = 0.3. 
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